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Top quark pair production and decay at hadron colliders: Predictions at 
NLO QCD including spin correlations* 

W. Bernreuther''^ A. Brandenburg''*, Z.G. Si""*, P. Uwer'= 

'^Institut f. Theoretische Physik, RWTH Aachen, 52056 Aachen, Germany 

'^DESY-Theorie, 22603 Hamburg, Germany 

■^Institut f. Theoretische Teilchenphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 

We present results at NLO QCD for hadronic production and decay of top quark pairs, taking into account tt 
spin correlations. 
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The interactions of the top quark have not been 
precisely tested so far. The expected data sam- 
ples of Run II of the Tevatron collider (eventu- 
ally about 10* it pairs per annum) and the LHC 
(~ 10*" it pairs p. a.) will make possible detailed 
investigations of top quark production and decay. 
In particular, one may search for possible devia- 
tions from the QCD production mechanism and 
the V-A structure of the top decay by analysing 
energy and angular distributions of the tt decay 
products. Such studies will be feasible since top 
quarks do not hadronize due to their very short 
lifetime, and perturbation theory gives a reliable 
description of top quark production and decay. 
To search for new physics effects, precise theo- 
retical predictions for top quark pair production 
and decay at hadron colliders within the Standard 
Model are needed. We consider here in particular 
the following processes: 



(1) 



where h\ i — p,p; I — e,^,T, and jt [ji) denote 
jets originating from hadronic t (t) decays. An 
observable that is very sensitive to the dynamics 




*Report DESY 02-144, PITHA 02/15, TTP02-24 
tsupported by BMBF contract 05 HTl PAA 4 

Speaker at the conference. Supported by a Heisenberg 
fellowship of D.F.G. 



of top quark production and decay in the above 
reactions is the double differential angular distri- 
bution of the top decay products, e.g. for the 
dilepton channel: 
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(1 - C cos 6*+ cos 6l_). (2) 



In (H) , 9± are the angles between the £^ direction 
of flight in the t{t) rest frame with respect to ar- 
bitrary axes a (b) which will be specified below. 
(Terms linear in cos9± are absent in (|^) for our 
choices of a, b due to parity invariance of QCD.) 

The calculation of the distribution (||) at next- 
to-leading order (NLO) QCD simplifies enor- 
mously in the leading pole approximation (LPA), 
which amounts to expanding the full amplitudes 
for (|^) around the complex poles of the t and 
i propagators. Only the leading pole terms are 
kept in this expansion, i.e. one neglects terms of 
order Tt/rrit « 1%. Within the LPA, the radia- 
tive corrections can be classified into factorizable 
and non-factorizable [Q contributions. Here we 
consider only the factorizable corrections. In this 
case the coefficient in the distribution (0) factor- 



izes: 

C = K_|_K__D. 

In (H), is the tt double spin asymmetry 
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(4) 
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where -/V(tT) denotes the number of tt pairs with 
t {t) spin parallel to the reference axis a (b) etc. 
The directions a and b can thus be identified with 
the spin quantization axes for t and t, and D re- 
flects the strength of the correlation between the 
t and t spins for the chosen axes. Further, k± is 
the spin analysing power of the charged lepton in 
the decay t{t) b(b)£^i'{D) defined by the decay 
distribution 



1 dr 



1± K± cos-d± 
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where d± is the angle between the t (t) spin and 
the direction of flight. From M we obtain 



1 - 0.015as 



(6) 



which means that the charged lepton serves as 
a perfect analyser of the top quark spin. For 
hadronic decays t bqq' one has a decay distri- 
bution analogous to (j^), and the spin analysing 
power of jets can be defined. To order as (^] (and 
using as{mt) — 0.108), 

Kb = -0.408 X (1 - O.UOas) = -0.393, (7) 
Kj = +0.510 X (1 - 0.654a,) = +0.474, (8) 

where Kj is the analysing power of the least en- 
ergetic non-b-quark jet. The results for C given 
below are for the dilepton channel. The corre- 
sponding results for the single lepton channel are 
obtained by replacing K4. (or k„) in Eq. (^) by ki, 
or Kj . The loss of analysing power using hadronic 
final states is overcompensated by the gain in 
statistics and in efhciency to reconstruct the t (t) 
rest frames. 

To compute D at NLO QCD, the differential 
cross sections for the following parton processes 
are needed to order a,, where the full information 
on the t and t spins has to be kept: 



qq 



tt, ttg: gg^tt, ttg; q{q)g ttq{q). (9) 



Results at NLO QCD for the MS subtracted par- 
ton cross sections a for the above processes with 
ti spins summed over have been obtained in 
§, while 
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Figure 1. Dimensionless scaling functions g^^\ri) 
(full), g^'^\r]) (dotted), and g^'^H'n) (dashed) that 
determine aD in the helicity basis for the gg ini- 
tial state. 

has been computed for different spin quantiza- 
tion axes in |Q. This combination of spin- 
dependent cross sections can be written as follows 

^,2 r 



aD = 



with 



g(0)(77)+4™,g« 



gW=gWi^)+g(^)irj)ln{I^ 



(11) 



(12) 



<7i? = a(tT) + 'T(ii) - <7(Ti) - adt) 



(10) 



where 77 = s/{4m'^) — 1, and we use a common 
scale /i for the renormalization and factorization 
scale. As an example, Fig. 1 shows, for the parton 
process gg — > ti{g), the three functions that de- 
termine aD for the helicity basis, i.e. a (b) is cho- 
sen to be the t (t) direction of flight. Apart from 
the helicity basis we also consider the beam basis, 
where a and b are equal to the direction of one of 
the hadron beams in the laboratory frame, and 
the so-called off-diagonal basis [||, which is de- 
fined by the requirement that a{1l) = <^{l^) = 
for the process qq tt at tree level. 

Table 1 lists our results ||] for the coefficient 
C in the double differential distribution (0) for 
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Table 1 

Coefficient C defined in Eq. (^) to leading and 
next-to-leading order in for the parton distri- 
bution functions of ||T^ and fi ~ mt ^ 175 GeV. 





PP, = 


2 TeV 


PP, Vs = 


14 TeV 




LO 


NLO 


LO 


NLO 


Choi. 


-0.456 


-0.389 


0.305 


0.311 


Cbeam 


0.910 


0.806 


-0.005 


-0.072 


Coff. 


0.918 


0.813 


-0.027 


-0.089 



Table 2 

Upper part: Dependence on /i of the correlation 
coefficients C, computed at NLO with the PDF 
ofref. Lower part: Correlation coefficients C 
at NLO for fi = rat and different sets of parton 
distribution functions: GRV98[|l), CTEQ5|0, 
and MRST98 (c-g)|l|. 







Tevatron 




LHC 




Choi. 


Cbcam 


Coff'. 


Choi. 


mt/2 


-0.364 


0.774 


0.779 


0.278 


mt 


-0.389 


0.806 


0.813 


0.311 


2mt 


-0.407 


0.829 


0.836 


0.331 


PDF 


Choi. 


C*bcain 


Coff. 


Chcl. 


GRV98 


-0.325 


0.734 


0.739 


0.332 


CTEQ5 


-0.389 


0.806 


0.813 


0.311 


MRST98 


-0.417 


0.838 


0.846 


0.315 



these different choices of spin quantization axes, 
using the parton distribution functions (PDFs) of 
10 1 . One sees that at the Tevatron, the spin cor- 
relations are largest in the beam and off-diagonal 
basis. The QCD corrections reduce the LO results 
by about 10%. For the LHC, the best choice of 
the above three bases is the helicity basis, and the 
QCD corrections are small in that case. 

Table 2 shows the dependence of the NLO re- 
sults for C on the scale ^ (upper part) and on the 
choice of the PDFs (lower part) . At the Tevatron 
the spread of results for different PDFs is larger 
than the scale uncertainty. This rather strong 
dependence on the PDFs is largely due to the 
fact that the contributions from qq and gg initial 
states contribute to C with opposite signs. This 
may offer the possibility to constrain the quark 



and gluon content of the proton by a precise mea- 
surement of the double angular distribution (||). 

In all results above we used mt = 175 GeV. A 
variation of nit from 170 to 180 GeV changes the 
results for the Tevatron, again for fi — mt and 
PDFs of ^ as follows: Chei. varies from -0.378 
to -0.397, Cbcam from 0.790 to 0.817, and CoS. 
from 0.797 to 0.822. For the LHC, Choi, changes 
by less than a percent. 

The results above have been obtained without 
imposing any kinematic cuts. Standard cuts on 
the top quark transverse momentum and rapidity 
only have a small effect on C: For the Tevatron, 
demanding |k^j| > 15 GeV and \rtj\ < 2 leads to 
the following results: Chei. = —0.386, Cbcam = 
0.815, Coff. = 0.823. For the LHC, when impos- 
ing the cuts |k^j| > 20 GeV and \rt^t\ < 3, we 
find Chei. = 0.295. 

In summary, tt spin correlations are large ef- 
fects that can be studied at the Tevatron and 
the LHC by measuring double angular distribu- 
tions both in the dilepton and single lepton decay 
channels. The QCD corrections to these distribu- 
tions have been computed and are under control. 
Spin correlations are suited to analyse in detail 
top quark interactions, search for new effects, and 
may help to constrain the parton content of the 
proton. 
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